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Abstract 


Results  are  presented  from  ab  initio  quantum  chemical  calculations  on  the  ground  state  potential 
energy  surface  for  the  open  shell  doublet  reaction: 

HNO  +  NO  -  N20  +  OH.  (I) 

This  reaction  is  thought  to  be  important  in  the  sustained  combustion  of  nitrogen-containing  fuels. 
Geometry  optimizations  have  been  carried  out  at  the  QCISD/6-31G(d,p)  level  using  UHF  zeroth- 
order  wave  functions  in  search  of  local  minima  and  transition  states,  followed  by  energy  refinements 
for  all  critical  points  at  the  QCISD(T)/6-31 1  +G(2df,2p)  level.  All  saddle  points  were  subjected  to 
intrinsic  reaction  coordinate  calculations  to  determine  the  minima  to  which  they  connect.  We  report 
two  completed  reaction  pathways  for  (I).  Some  of  these  results  will  be  compared  with  earlier  work 
by  others  on  the  same  reaction  where  geometry  optimizations  were  carried  out  primarily  at  the 
UMP2/6-311G(d,p)  level.  The  comparison  brings  to  light  some  interesting  differences  between 
QCISD  and  UMP2  theory  in  the  prediction  of  molecular  structures  and  the  connectivity  of  critical 
points  (and  therefore,  connected  reaction  paths)  for  the  doublet  system  in  Reaction  I. 
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1.  Introduction 


Chemical  reactions  involving  HNO  and  NO  are  considered  to  be  key  processes  in  the  combustion 
of  nitrogen-containing  fuels.  Our  primary  interest  in  NO  chemistry  centers  upon  its  role  in 
controlling  ignition  chemistry  in  the  dark  zone  of  nitramine-based  gun  propellants.  One  of  the 
distinctive  features  observed  during  the  combustion  of  nitramine-containing  solid  propellants  is  the 
two-stage  nature  of  the  flame.  At  pressures  of  ~  2  MPa,  there  is  a  separation  between  the  solid 
surface  and  the  gas  phase  luminous  flame  that  is  typically  0.5  cm  wide.  This  nonluminous  region 
is  commonly  referred  to  as  the  dark  zone.  The  dark  zone  is  believed  to  contain  NO  and  NO- 
containing  intermediates  that  react  slowly.  Detailed  kinetic  modeling  of  the  dark  zone  coupled  with 
a  sensitivity  analysis  of  the  reaction  set  used  in  the  model  by  Vanderhoff  and  coworkers  [  1  ]  indicate 
that  the  reaction 


HNO  +  NO  -  N20  +  OH  (I) 

is  the  most  sensitive  in  controlling  the  ignition  chemistry  of  nitramine  propellants.  This  reaction  is 
the  rate-controlling  step  in  the  thermal  reduction  of  NO  by  H2  at  temperatures  greater  than  1,000  K. 
In  addition,  this  reaction  is  also  known  to  initiate  the  formation  of  N2and  H20  [2, 3]. 

There  have  been  theoretical  studies  of  related  reactions  involving  hydrogen  and  NO,  where  some 
of  the  more  recent  work  includes  the  reactions  HNO  +  H  [4, 5],  HNO  +  OH  [6],  and  NH  +  N02  [7]. 
To  date,  there  has  been  no  direct  experimental  measurement  of  the  rate  for  the  HNO  +  NO  reaction, 
presumably  due  to  the  difficulties  associated  with  making  HNO.  The  rate  constant,  however,  has 
been  estimated  by  Wilde  [8]  to  be  10 123 1 03  exp(-260,000  ±  5,000/RT)  cm3/(mol-s)  over  the 
temperature  range  800-1 ,060°  C.  This  estimate  was  obtained  from  kinetic  modeling  by  Graven  [9], 
and  from  measurements  for  the  rate  of  NO  disappearance  by  Kaufman  and  Decker  [10]  and  pressure 
changes  by  Hinshelwood  et  al.  [1 1,  12],  all  of  whom  were  studying  the  kinetics  of  the  H2  +  NO 
system. 
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In  light  of  the  suggestion  that  Reaction  I  is  perhaps  the  rate-controlling  step  in  the  chemistry  of 
the  dark  zone  of  a  nitramine  propellant  flame,  we  initiated  an  ab  initio  study  of  the  reaction  of  HNO 
+  NO.  During  the  preliminary  stages  of  our  calculations,  Mebel,  Morokuma,  Lin,  and  Melius 
(MMLM)  [13]  published  many  critical  points  for  the  HNO  +  NO  potential  energy  surface  (PES)  and 
information  about  the  reaction  paths  leading  to  products.  All  of  these  critical  points  were  described 
at  the  UQCISD(T)/6-31  lG(d,p)//UMP2/6-31  lG(d,p)  level.  Energies  of  the  critical  points  on  their 
“minimum  energy  reaction  path  (MERP)”  [13]  were  refined  using  a  modified  version  of  the 
Gaussian-2  (G2)  [14-16]  scheme  [referred  to  as  G2(PU)].  In  their  G2(PU)  approach, 
UMP2/6-311G(d,p)  [17]  (i.e.,  projected  UHF-MP2  [17])  is  used  to  optimize  all  geometries  (with 
two  exceptions).  The  geometries  of  reactants  and  the  transition  state  corresponding  to  the 
rate-determining  step  were  optimized  at  the  UQCISD/6-31  lG(d,p)  level  and  used  in  the  G2(PU) 
scheme.  The  results  of  their  study  suggest  that  the  rate-determining  barrier  for  Reaction  I  is  a  1,3 
hydrogen  migration  with  a  barrier  height  of  21.6  kcal/mol.  In  an  earlier  study,  Mebel,  Morokuma, 
and  Lin  [7]  used  UMP2  theory  to  predict  structures  on  the  PES  of  a  related  reaction  where 
NH  +  N02  proceeds  to  two  products,  namely,  HNO  +  NO  and  N20  +  OH,  which  are  the  reactants 
and  products  of  Reaction  I,  respectively. 

Our  early  results  on  this  system  led  us  to  believe  that  UMP2  is  not  capable  of  adequately 
describing  the  electronic  structure  in  this  system.  This  conclusion  is  based  on  the  fact  that  many 
points  on  the  PES  of  Reaction  I  have  UHF  zeroth-order  wave  functions  with  large  spin 
contamination,  where  some  S2  values  are  as  large  as  1 .36  (the  correct  value  is  0.75).  It  is  well  known 
that  open  shell  molecules  described  by  UHF  orbitals  are  subject  to  such  spin  contamination.  This 
has  been  cited  as  one  possible  cause  for  the  slow  convergence  in  the  MP  series  [  1 8-2 1  ] .  In  choosing 
a  more  suitable  (yet  computationally  feasible)  approach,  we  decided  to  locate  critical  points  on 
portions  of  the  PES  at  the  QCISD  level.  Studies  have  shown  that  coupled  cluster  theory,  of  which 
QCISD  is  a  subset,  is  effective  in  correcting  for  spin  contamination  [22]  in  the  UHF  wave  function, 
and  QCISD  has  been  shown  to  predict  relative  energies  in  good  agreement  with  experiment  for 
another  doublet  system  [23].  As  further  justification  for  the  use  of  these  methods  in  predicting 
relative  energies  on  an  open  shell  system,  we  cite  the  work  of  Lee,  Rendell,  and  Taylor  [24]  on  Be 
and  Mg  clusters.  Lee,  Rendell,  and  Taylor  [24]  have  shown  that  relative  energies,  in  the  form  of 
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binding  energies,  calculated  at  the  CCSD(T)  and  QCISD(T)  levels  for  small  Be  and  Mg  clusters 
reproduce  .  .  the  accurate  MRCI  values  very  well.”  These  are  clusters  with  electronic  states 
purported  to  be  of  multireference  character  [24].  This  implies  that  the  approximate  inclusion  of 
triples  in  CCSD(T)  and  QCISD(T)  provides  some  of  the  important  correlation  corrections  obtained 
in  a  multireference  treatment.  For  the  above  reasons,  all  critical  points  on  the  PES  have  been 
optimized  at  the  UQCISD/6-31G(d,p)  level,  followed  by  energy  refinement  at  the 
UQCISD(T)/6-311+G(2df,2p)  level.  Basis  set  dependence  has  been  investigated  by  optimizing 
several  critical  points  on  the  PES  at  the  UQCISD  level  using  a  triple-zeta  basis  set. 

Most  of  the  critical  points  described  in  this  work  are  features  of  two  completed  reaction  paths 
for  Reaction  I.  These  paths  were  established  through  intrinsic  reaction  coordinate  (IRC)  calculations 
leading  from  the  various  transition  states  to  local  minima.  While  we  make  no  claim  to  having  found 
all  possible  pathways  for  Reaction  I,  the  results  presented  here  offer  a  convincing  argument  that  the 
two  paths  determined  in  this  study  are  continuously  connected  from  reactants  to  products.  The 
minimum  energy  pathway  reported  here  differs  from  that  reported  by  MMLM.  These  differences 
in  features  along  the  reaction  paths  could  lead  to  different  conclusions  regarding  mechanisms  for 
Reaction  I.  Significant  inconsistencies  are  noted  between  a  number  of  the  UQCISD  structures  from 
this  study  and  the  UMP2  structures  from  previous  works  [7, 13]. 

2.  Methods 

All  calculations  describing  the  PES  were  performed  using  the  Gaussian  92  [25]  (G92)  or 
Gaussian  94  [26]  (G94)  set  of  computer  codes.  Critical  points  on  the  PES  were  determined  through 
UQCISD  [27— 29]/6-3 1  G(d,p)  [30-3 1  ]  geometry  optimizations  (using  a  frozen  core),  where  gradients 
were  converged  to  at  least  the  default  settings.  Henceforth,  it  will  be  understood  that  all  zeroth-order 
wave  functions  are  of  the  UHF  type,  and  the  “U”  label  will  be  dropped  from  the  correlation  treatment 
designation.  Each  critical  point  was  characterized  through  normal  mode  analysis.  In  addition,  a 
stability  check  [32]  (with  default  perturbations)  was  run  on  the  UHF  wave  functions.  All  transition 
states  reported  here  had  one  imaginary  frequency,  and  all  minima  had  no  imaginary  frequencies.  The 


3 


transition  states  were  subjected  to  IRC  [33-34]  calculations  (using  the  default  step  size)  to  facilitate 
connection  with  minima  along  the  reaction  paths.  Each  IRC  terminated  upon  reaching  a  minimum 
as  defined  by  the  default  criteria  provided  in  G92  or  G94,  unless  otherwise  stated.  For  those  points 
that  did  not  terminate  in  this  fashion,  the  last  converged  point  on  the  reaction  path  was  used  as  the 
initial  structure  for  a  full  geometry  optimization.  For  each  instance  in  which  such  an  optimization 
was  attempted,  a  local  minimum  was  found.  It  is  assumed  that  this  is  the  minimum  to  which  the  IRC 
would  have  walked  if  convergence  problems  had  not  been  encountered.  Energies  and  geometries 
will  be  given  to  justify  the  conclusion  that  the  IRC  was  proceeding  to  the  same  minimum  located  by 
the  subsequent  optimization.  Once  all  critical  points  for  the  reaction  paths  were  determined,  the 
energies  were  refined  (using  a  frozen  core)  at  the  QCISD(T)  [27]/6-3 1 1  +G(2df,2p)  [35-38]//QCISD/ 
6-31G(d,p)  level.  These  refined  energies  were  then  corrected  for  zero-point  energies  (ZPE)  using 
the  QCISD/6-31G(d,p)  structures  and  frequencies. 

Some  critical  points  on  the  PES  have  spin  S2  values  substantially  above  the  0.75  magnitude 
expected  for  a  doublet  system.  The  proper  way  to  assess  the  capability  of  the  QCISD  method  for 
correcting  this  spin  contamination  in  the  UHF  wave  function  is  to  calculate  the  S2  expectation  value 
over  the  QCISD  density.  Unfortunately,  we  do  not  have  access  to  codes  that  include  this  capability. 
We  can,  however,  calculate  this  expectation  value  for  the  coupled-cluster  singles  and  doubles 
[39-41]  (CCSD)  method  (with  no  frozen  orbitals)  using  our  version  of  the  ACES  II  quantum 
chemistry  package  [42-43],  which  includes  the  code  for  <S2>CC  written  by  Stanton  [22].  Even 
though  QCISD  is  only  a  subset  of  the  full  representation  of  the  CCSD  equations  [24, 35, 36, 44-46] , 
it  will  be  assumed  that  the  two  approaches  give  similar  expectation  values  of  S2. 

3.  Results  and  Discussion 


Absolute,  relative,  and  ZPE  and  S2  values  of  the  critical  points  on  the  HNO  +  NO  PES  are  given 
in  Table  1.  Harmonic  vibrational  frequencies  of  all  critical  points  calculated  at  the  QCISD/6- 
3  lG(d,p)  level  are  given  in  the  appendix  (Table  A-l).  We  have  calculated  the  CCSD  S2  expectation 
values  for  five  selected  critical  points  on  the  PES  where  the  UHF  S2  is  large.  From 
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Table  1.  Absolute  Energies,  ZPE,  Zero-Point  Corrected  Relative  Energies,  and  S2  Values  for  Species  on  the  HNO  +  NO  PES.  All 
Geometries  Optimized  at  QCISD/6-31G(d,p). 
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Table  1  we  see  that  the  S2,s  before  and  after  the  CCSD  calculations  are  A(1.36  -  0.81); 
E(1.28  -*  0.77);  F(1.36  -  0.82);  L(1.28  -*  0.77);  N(l.ll  -*  0.76).  The  most  heavily  contaminated 
UHF  wave  function  has  its  spin  corrected  to  0.82,  where  the  correct  value  is  0.75  for  a  doublet 
system. 

It  has  been  shown  that  QCISD  neglects  certain  terms  in  the  energy  expression  produced  by  the 
single  excitation  operator  Tx  that  are  present  in  the  complete  CCSD  treatment  [24, 35, 36,  44-46]. 
Therefore,  QCISD  could  produce  significantly  different  relative  energies  along  the  reaction  paths 
than  CCSD,  which  would  cast  doubt  upon  the  ability  of  QCISD  to  achieve  the  primary  goal  of  this 
study,  which  is  to  accurately  map  the  relative  energetics  defining  the  PES.  Lee,  Rendell,  and  Taylor 
[24]  concluded  from  their  study  of  18  different  chemical  species  that  “. . .  the  better  the  SCF  wave 
function  approximates  the  true  wave  function,  the  better  the  agreement  between  the  QCISD  and 
CCSD  correlation  energies.”  To  ascertain  the  impact  on  the  relative  energies  due  to  neglecting  the 
Trcontaining  terms,  CCSD  [42-43]/6-31G(d,p)  calculations  were  performed  on  the  five  points 
mentioned  above  with  the  worst  UHF  S2  values.  Table  2  contains  the  energies  of  each  of  the  five 
structures  relative  to  the  reactants  for  QCISD  vs.  CCSD.  This  table  also  contains  the  absolute  total 
energies  at  each  point  for  both  QCISD  and  CCSD  methods  and  the  difference  between  the  absolute 
QCISD  and  CCSD  energies  at  each  point.  All  calculations  were  done  with  frozen  core  orbitals  to 
be  consistent  with  our  geometry  optimizations.  The  geometries  and  ZPE  corrections  were  taken 
from  the  QCISD/6-3 1  G(d,p)  calculations.  There  are  several  points  worth  noting  concerning  Table  2. 
As  pointed  out  by  Lee  et  al.  [24],  QCISD  tends  to  overestimate  the  correlation  correction  and  hence 
produces  absolute  energies  lower  than  those  of  CCSD.  This  trend  is  again  seen  in  the  data  of 
Table  2.  The  first  important  point  to  note  is  that  the  difference  between  the  QCISD  and  CCSD 
absolute  energies  for  each  structure  “i  ”  defined  as  AE^5  =  E^QCISD)  -  E,(CCSD),  is  nearly 
constant  over  the  five  points  (A,  E,  F,  L,  N)  on  the  PES,  falling  in  the  range  0.00484-0.00528 
hartrees  or  3.04-3.31  kcal/mol.  The  second  important  point  is  that  the  difference  in  QCISD  vs. 
CCSD  energies  for  the  reactants  is  2.59  kcal/mol.  For  the  five  structures  on  the  PES,  this  produces 
a  maximum  and  average  difference  in  relative  energies  between  CCSD  and  QCISD  of  0.72  and 
0.56  kcal/mol,  respectively.  The  largest  difference  belongs  to  structure  F,  which  also  has  the  worst 
UHF  S2  value.  Based  on  this  data,  it  is  assumed  that  the  QCISD  PES  for  this  reaction  would  agree 
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Table  2.  QCISD76-31G(d,p)  vs.  CCSD/6-31G(d,p)  Energies  for  Selected  Structures  From  the  PES  With  Large  S2  Values. 
QCISD/6-31G(d,p)  Geometries  and  ZPE  Corrections  Used  for  All  Structures. 
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with  a  complete  CCSD  PES  to  within  ~1.0  kcal/mol  or  less  for  all  structures,  which  is  probably 
smaller  than  the  error  in  the  AE’s  associated  with  this  level  of  theory.  Hence,  no  attempt  was  made 
to  repeat  these  computationally  intensive  calculations  at  the  CCSD  level  for  the  remaining  points  on 
the  PES. 

T wo  paths  leading  to  formation  of  the  products  of  Reaction  I  have  been  established  through  IRC 
calculations  from  the  transition  states  optimized  at  the  QCISD/6-31G(d,p)  level.  Figures  1  and  2 
show  diagrams  of  relative  energies  along  reaction  Paths  1  and  2,  respectively.  The  geometries  of  the 
stationary  points  are  placed  at  their  appropriate  locations  along  the  reaction  paths.  Each  structure 
is  labeled  alphabetically  in  these  figures  to  associate  the  structure  with  the  species  listed  in  Tables  1 
and  2  and  in  the  following  discussion.  Unless  otherwise  stated,  all  barrier  heights  discussed 
throughout  this  report  correspond  to  the  highest  level  of  calculation  (QCISD(T)/6-3 1 1  +G(2df,2p)// 
QCISD/6-31G(d,f))  and  have  been  corrected  for  ZPE. 

The  two  reaction  pathways  are  similar  in  that  the  first  step  involves  formation  of  the  HN(0)N0 
intermediate  through  N-N  bond  formation,  which  then  undergoes  a  hydrogen  migration  as  the 
reaction  proceeds.  All  references  to  cis-  or  trans-  conformations  are  relative  to  the  four  heavy  atoms 
unless  otherwise  noted.  Path  1  (Figure  1),  whose  first  step  is  formation  of  cis  -HN(0)NO,  proceeds 
to  products  through  a  1 ,2  hydrogen  migration  reaction.  Path  2  (Figure  2)  proceeds  through  formation 
of  trans-HN(0)NO,  which  subsequently  undergoes  a  1 ,3  hydrogen  shift.  Our  results  begin  to  differ 
from  MMLMs  immediately  after  the  formation  of  the  entrance  channel  TS.  In  our  study,  the 
cis-entrance  channel  barrier  A  in  Path  1  goes  on  to  form  a  minimum  cis-structure  B.  The  trans¬ 
entrance  channel  TS  F  in  Path  2  remains  trans-  as  it  proceeds  to  the  minimum  structure  G.  MMLM 
report  a  single  cis-entrance  channel  transition  state,  1 ',  that  proceeds  to  a  minimum  trans-structure,  2. 
For  comparison,  the  MERP  of  MMLM  has  been  reconstructed  in  Figure  3,  with  the  reported 
structures  (including  1 '  and  2)  at  each  critical  point.  At  this  point,  we  will  discuss  features  of 
reaction  Paths  1  and  2  separately. 

3.1  Path  1.  Critical  point  A  is  a  loose,  nonplanar  transition  state  with  an  N-N  distance  of  1 .90  A 
and  HNO  and  NO  moieties  closely  resembling  the  isolated  molecules  (see  Figure  1).  The 
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Pathway  Connects  to  Pathway  2. 


Figure  2.  Pathway  2.  Reported  Here  as  the  Minimum  Energy  Reaction  Pathway  for  Reaction  I.  AI1  Structures  and  Relative 
Energies  Are  for  Zero-Point  Corrected  Energies  Calculated  at  the  QCISD/6-311+G(2df,2p)//QCISD/6-31G(d,p)  Level. 


intermediate  formed  upon  crossing  this  barrier  (structure  B  in  Figure  1)  is  cis-planar  and  has  an  N-N 
bond  distance  of  1.42  A.  IRC  calculations  for  the  walk  A  -  B  terminated  upon  reaching  the  local 
minimum  B.  The  IRC  walk  of  A  -  REACTANTS  terminated  when  the  energy  was  within 
0.43  kcal/mol  of  the  reactant  asymptote  and  the  geometric  parameters  of  the  species  were  within 
0.2%  or  less  of  the  reactants’  values.  Two  reaction  pathways  were  found  starting  from  B.  The  lower 
energy  path  proceeds  over  the  TS  structure  representing  a  simple  torsion  around  the  N-N  bond.  IRCs 
from  this  TS,  denoted  as  structure  M  in  Figure  1,  connect  Paths  1  and  2  by  the  relationship 
B  -  M  -  G. 

Barrier  M  to  isomerization  from  B  to  G  is  14.1  kcal/mol  relative  to  B.  The  nitrogen  to  which 
hydrogen  is  bonded  is  nonplanar  with  a  torsion  angle  of  t(N4-N2-03-H1)  =  122°.  The  features  of 
this  portion  of  the  path  suggest  that  the  reactants  might  initially  start  out  on  one  reaction  path  and 
move  easily  to  the  other  through  the  isomerization  from  B  to  G  (or  vice  versa)  to  continue  on  the 
other  reaction  path.  Based  upon  the  comparative  heights  of  barriers  M  (13.4  kcal/mol)  and  F  (the 
entrance  channel  barrier  on  Path  2)  (1 1.3  kcal/mol)  relative  to  reactants,  the  formation  of  structure 
G  is  slightly  favored  energetically  through  traversing  the  entrance  barrier  F.  However,  under 
combustion  conditions,  the  isomerization  barrier  M  is  also  accessible,  so  Path  2  could  be  reached 
by  either  TS  A  (via  M)  or  F.  Neither  MMLM  nor  Mebel,  Morokuma,  and  Lin  [7]  report  this 
transition  state  M. 

The  TS  most  closely  resembling  M  at  the  UMP2/6-311G(d,p)  level  was  given  by  Mebel, 
Morokuma,  and  Lin  [7]  as  a  planar  structure,  labeled  6'.  The  planarity  of  6'  is  in  sharp  contrast  with 
the  QCISD  prediction  of  109°  for  the  torsion  angle  x(ONNO).  An  attempt  was  made  to  locate  this 
planar  transition  state  6'  using  the  basis  set  in  the  Mebel,  Morokuma,  and  Lin  [7]  study 
[i.e.,  6-31  lG(d,p)]  but  at  the  QCISD  level  of  correlation.  Before  the  QCISD/6-31  lG(d,p)  TS  search 
was  attempted,  we  reproduced  planar  structure  6'  from  Mebel,  Morokuma,  and  Lin  [7]  through  a 
transition  state  optimization  and  subsequent  normal  mode  analysis  at  the  MP2/6-31  lG(d,p)  level. 
The  MP2  geometry  (which  agreed  with  that  published  [7])  and  corresponding  force  constant  matrix 
(FCM)  were  used  as  the  initial  conditions  for  the  QCISD/6-31  lG(d,p)  TS  search.  The  QCISD 
transition  state  search  did  NOT  locate  a  cis-planar  saddle  point,  rather  it  found  a  cis-planar  minimum 
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that  is  essentially  identical  to  our  minimum  structure  B.  In  Figure  4,  the  QCISD/6-31  lG(d,p) 
structures  are  compared  with  the  purported  MP2/6-31  lG(d,p)  analogues  from  MMLM  and  Mebel, 
Morokuma,  and  Lin  [7].  Table  3  contains  the  relative  energies  for  these  structures,  and  Table  A-2 
reports  their  vibrational  frequencies.  In  spite  of  the  larger  basis  set,  all  geometrical  parameters  from 
the  QCISD/6-31  lG(d,p)  calculations  agree  with  the  QCISD/6-31G(d,p)  to  within  1%  for  B.  As  a 
final  attempt  to  locate  6 '  at  the  QCISD/ 6-3 1  lG(d,p),  the  FCM  was  first  calculated  using  the  reported 
structure  in  Mebel,  Morokuma,  and  Lin  [7].  This  FCM  and  the  Mebel,  Morokuma,  and  Lin  [7] 
structure  were  then  used  as  the  initial  guess  in  a  TS  search,  which  again  converged  on  our  local 
minimum  B.  We  were  unable  to  locate  the  Mebel,  Morokuma,  and  Lin  [7]  TS  6'  at  the  QCISD  level. 
The  near  exact  agreement  between  the  QCISD  predictions  using  the  6-31G(d,p)  and  6-311G(d,p) 
basis  sets  suggests  that  the  use  of  the  smaller  double-zeta  basis  adequately  describes  the  valence 
space  with  respect  to  predicting  reliable  molecular  structures  for  this  system,  at  least  for  the 
minimum  structures. 

Returning  to  Path  1  (shown  in  Figure  1),  the  reaction  proceeds  from  B  to  D  through  a  TS 
representing  a  1,2  hydrogen  shift  to  03.  This  transition  state  is  labeled  C  in  Figure  1.  The  barrier 
to  this  hydrogen  migration  (5 1 .2  kcal/mol)  is  the  largest  along  this  path.  Therefore,  this  1 ,2  H-shift 
is  the  rate-determining  step  on  Path  1 .  The  magnitude  of  this  barrier  precludes  this  reaction  path  as 
a  contributor  to  thermal  reaction  under  combustion  conditions.  While  many  steps  were  completed 
in  the  IRCs  going  from  C  toward  B  and  D,  none  went  to  completion.  The  number  of  successful  steps 
for  each  IRC  is  C  -  B(18);  C  -  D(14).  Optimizations  were  then  performed  using  the  geometries 
corresponding  to  the  last  successful  IRC  step  in  each  walk,  and  the  resulting  structures  were 
consistent  with  B  and  D,  respectively. 

MMLM  report  a  planar  saddle  point  structure  resembling  our  C  (labeled  4"  in  their  paper),  which 
also  connects  minima  similar  to  B  and  D  through  a  1 ,2  hydrogen  shift.  They  do  not,  however,  report 
any  nonplanar  structure  akin  to  our  C.  A  TS  search  was  initiated  at  the  QCISD/6-31  G(d,p)  level  to 
locate  a  planar  analogue  of  C.  Such  a  structure  was  located  and  is  labeled  C'  in  Figure  1. 
Vibrational  analysis  supports  the  assignment  of  this  structure  as  a  TS.  Subsequent  IRC  calculations 
did  not  complete  successfully,  but  ensuing  optimizations  connect  this  TS  to  B  and  D.  The  number 
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Table  3.  Geometry  Optimizations  at  QCISI)/6-311G(d,p):  Absolute  Energies,  Zero-Point  Energies  (ZPE) ,  Zero-Point  Corrected 
Relative  Energies,  and  S2  Values  for  Selected  Species  on  the  HNO  +  NO  PES 
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of  successfully  converged  steps  along  each  IRC  prior  to  optimization  are  C'  -  B(17);  C'  -  D(18). 
The  energy  of  C  and  C'  relative  to  B  are  nearly  identical,  with  C'  lower  than  C  by  less  than 
1  kcal/mol. 

To  eliminate  the  possibility  that  our  nonplanar  C  exists  only  as  an  artifact  of  the  smaller  basis 
set  used  in  this  study,  C  was  again  sought  at  the  QCISD/6-311G(d,p)  level  starting  with  the 
QCISD/6-3  lG(d,p)  optimized  structure  and  force  constant  matrix.  This  converged  on  a  nonplanar 
TS  verified  through  vibrational  analysis.  The  resulting  structure  had  bond  lengths  that  differed  by 
1  %  or  less  from  the  structure  calculated  with  the  double-zeta  basis  set,  the  simple  angles  differed  by 
less  than  0.6%,  and  the  largest  torsional  difference  was  only  3  ° .  Once  again,  the  6-3 1  G(d,p)  and  6- 
3 1  lG(d,p)  calculations  gave  quantitatively  similar  structures  at  the  QCISD  level,  this  time  for  a  TS. 

Finally,  a  saddle  point  search  was  conducted  to  locate  the  planar  C'  structure  at  the 
QCISD/6-3 1  lG(d,p)  level.  This  search  converged,  but  the  resulting  vibrational  analysis  produced 
two  imaginary  eigenvalues  (i.e.,  2057/  and  173/  cm'1).  Repeated  attempts  to  locate  a  true  TS  for 
the  planar  structure  with  the  triple-zeta  basis  failed.  These  results  cast  doubt  on  the  existence  of  the 
planar  TS  C'.  It  is  also  possible  that  a  1,3  hydrogen  shift  to  05  could  occur;  however,  no  attempt 
was  made  to  locate  a  saddle  point  for  this  type  of  hydrogen  migration. 

Once  barrier  C  (or  C')  is  traversed,  another  intermediate  is  formed  (structure  D,  Figure  1)  in 
which  the  H-0  bond  distance  is  close  to  the  asymptotic  product  value,  the  N-H  separation  is  large 
at  1.87  A,  and  the  N-0  is  elongated  to  1.42  A.  A  final  saddle  point,  E,  leading  to  products  has  a 
structure  in  which  both  N-H  and  N-0  distances  are  long  (2.04  and  1.69  A,  respectively),  indicative 
of  the  O-H  separation  from  the  N20  moiety.  Also,  the  NNO  angle  is  approaching  that  of  the  linear 
conformation  for  the  separated  N20  product.  The  magnitude  of  this  barrier  relative  to  the  minimum 
for  structure  D  is  4.6  kcal/mol.  IRC  calculations  for  the  walk  E  -  D  proceeded  until  the  local 
minimum  corresponding  to  D  was  reached.  The  large  energy  release  upon  crossing  barrier  C 
(59  kcal/mol),  and  the  small  exit  channel  barrier  E  suggests  that  the  intermediate  D  would  be 
extremely  short-lived.  The  IRC  calculation  for  E  -  PRODUCTS  was  terminated  when  the  energy 
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was  within  0.28  kcal/mol  of  the  asymptote  and  the  geometric  parameters  were  within  0.3%  or  less 
of  products. 

3.2  Path  2.  Path  2  is  the  MERP  in  this  study,  which  can  be  compared  with  the  MERP  of 
MMT.M  as  shown  in  Figure  3.  The  entrance  channel  transition  state  structure,  F  (see  Figure  2), 
leading  to  formation  of  the  HN(0)N0  intermediate,  is  very  similar  to  the  entrance  channel  transition 
state  for  Path  1  except  for  a  trans-structure  in  the  heavy  atoms.  It  is  also  a  very  loose  TS  with  a  N-N 
separation  of  1.86  A.  The  barrier  height  is  11.3  kcal/mol,  which  is  2.9  kcal/mol  higher  than  the 
analogous  entrance  channel  barrier  for  Path  1.  The  IRC  walk  toward  reactants  was  halted  after  46 
steps.  The  structure  at  this  point  on  the  reaction  path  had  an  N-N  separation  of  2.84  A,  and  the 
energy  was  within  0.4  kcal/mol  of  the  isolated  reactants.  Also,  all  remaining  structural  parameters 
were  within  0.2%  of  the  reactants.  Upon  crossing  this  barrier,  a  minimum  trans-planar  HN(0)N0 
intermediate,  G,  is  formed  that  is  similar  to  its  cis-planar  counterpart  in  Path  1 .  The  IRC  walk  F  -  G 
aborted  after  nine  steps  due  to  the  updated  Hessian  having  the  “wrong  number  of  negative 
eigenvalues.”  All  of  the  bond  lengths  and  angles  agreed  to  within  0.2%  of  those  in  G,  and  the  largest 
deviation  between  the  two  structures  was  a  3°  difference  in  the  ONNO  dihedral  angle.  Also,  the 
energy  difference  between  this  last  point  on  the  path  and  G  was  only  0.03  kcal/mol.  Structure  G  is 
0.3  kcal/mol  lower  in  energy  than  its  cis-counterpart  B.  The  reader  is  reminded  that  structure  G 
connects  to  structure  B  on  Path  1  by  the  transition  state  M,  as  discussed  in  the  preceding  section. 

The  path  leading  from  F  to  G  represents  the  first  of  several  differences  in  the  connections  along 
the  MERP  of  this  study  and  that  of  MMLM.  MMLM  report  the  entrance  channel  transition  state  to 
their  MERP  to  be  cis-(structure  1 '),  which  then  walks  to  a  trans-minimum,  structure  2  (see  Figure  3). 
The  current  study  predicts  a  trans-entrance  channel  TS  (F)  proceeding  to  a  trans-minimum  structure 
G  (see  Figure  2).  The  nonplanar  structure  F  of  this  study  most  closely  resembles  the  MMLM  planar 
TS  structure  1.  For  comparison,  MMLM’s  structure  1  is  reproduced  in  Figure  4  along  with  our 
structure  F.  The  large  N-N  distance  strongly  suggests  a  very  weak  potential  for  rotation  about  the 
N-N  bond,  and  hence,  no  attempt  was  made  to  explore  the  issue  of  planarity  in  this  conformer. 
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Figure  3.  Reconstruction  of  MMLM’s  [13]  Minimum  Energy  Reaction  Pathway  for  Comparison  With  the  Current  Study’s 
Pathway  2.  Energies  Taken  From  Figure  2  of  Mebei  et  al.  [13]  for  Values  Along  the  Solid  Bold  Curve. 
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After  crossing  the  entrance  channel  barrier.  Path  2  now  proceeds  through  a  1,3  hydrogen  shift 
to  05.  The  barrier  (structure  H,  Figure  2)  for  the  1,3  shift  is  24.0  kcal/mol  relative  to  the  reactants 
(25.0  kcal/mol  relative  to  structure  G).  This  barrier  represents  the  rate-limiting  step  for  this  path, 
and  its  barrier  height  is  close  to  the  best  predictions  ofMMLM  (21.6  kcal/mol).  Wilde  [8]  extracted 
an  experimental  estimate  of  26  ±  5  kcal/mol  for  this  barrier  height  from  the  data  of  Kaufman  and 
Decker  [10]  and  Graven  [9],  above  1,000°  C.  This  experimental  estimate  is  in  good  agreement  with 
our  value  of  24.0  kcal/mol  for  H.  The  IRC  for  TS  H  toward  G  was  terminated  after  37  steps  along 
the  path.  At  this  point,  the  energy  was  within  0.01  kcal/mol  of  that  corresponding  to  structure  G. 
All  bond  lengths  and  angles  differed  by  no  more  than  0.2%  from  the  optimized  G  structure,  and 
torsion  angles  differed  by  less  than  1°.  The  IRC  from  TS  H  toward  products  led  to  a  local  energy 
minimum,  denoted  as  structure  I.  The  IRC  had  proceeded  to  a  structure  whose  energy  was  within 
0.04  kcal/mol  of  the  optimized  structure  I,  and  all  structural  parameters  differed  from  those  of 
structure  I  by  less  than  1%.  I  shows  formation  of  the  O-H  bond  concurrent  with  the  breaking  of  the 
N-H  bond.  The  H-0  bond  distance  is  close  to  the  asymptotic  product  value,  and  the  H-N  separation 
is  large  at  2.04  A. 

Structure  I  resides  in  a  very  shallow  well,  from  which  the  reaction  proceeds  over  barrier  J,  which 
is  1 .5  kcal/mol  higher  in  energy.  The  most  significant  change  in  geometry  as  TS  J  is  traversed  is  the 
opening  of  angle  N4-N2-03  as  it  approaches  the  linear  conformation  in  the  products.  The 
N4-N2-03  angle  of  the  final  intermediate,  K,  on  Path  2  has  opened  up  to  127°.  The  final  barrier  L 
represents  a  simultaneous  separation  of  the  O-H  group  from  the  N20  moiety  and  the  further  opening 
oftheN20  bond  angle  to  155°  as  it  approaches  the  linear  product.  This  final  barrier  is  11.3  kcal/mol 
relative  to  the  minimum  for  structure  K.  The  exit  channel  TS  along  MMLM’s  MERP  is  13.  TS  13 
(see  Figure  3)  corresponds  to  our  TS  structure  E  on  Path  1 .  MMLM  report  a  second,  higher  energy 
exit  channel  TS,  labeled  12  (not  shown  in  Figure  3),  which  has  a  cis-ONNO  arrangement,  while  TS 
L  is  trans.  It  is  certainly  possible  that  both  12  and  L  exist,  and  one  connects  to  the  other  by  a  simple 
rotation  about  the  N-N  bond.  The  barrier  to  rotation  is  likely  to  be  very  small,  and  the  molecule 
could  easily  pass  from  one  structure  to  the  other,  so  no  attempt  was  made  to  locate  TS  12. 
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3.3  Further  Comparisons  With  UMP2  Structures. 


3.3.1  J  •*  K.  Transition  state  structure  J  in  the  current  study  corresponds  most  closely  to 
MMLM’s  structure  14a  in  Figure  3.  The  torsion  angle  t(H1-05-N4-N2)  for  J  is  cis,  and  the 
QCISD/6-3 1  G(d,p)  IRC  walks  from  J  to  structure  K,  which  is  also  cis.  However,  in  MMLM’ s  study, 
the  IRC  from  this  transition  state  (structure  14a)  is  reported  to  walk  to  a  trans  (H-O-N-N)  planar 
structure  (7b,  Figure  3).  To  further  explore  this  contrast  between  the  two  studies,  we  attempted  to 
duplicate  the  calculations  of  MMLM  by  using  the  same  correlation  treatment  (MP2)  and  basis  set 
(6-3 1  lG(d,p)).  Using  the  default  tolerances  for  convergence  criteria  in  the  G92  set  of  programs,  and 
beginning  with  MMLM’s  published  structure  and  an  MP2  FCM,  we  located  the  same  transition  state 
structure  labeled  14a.  A  frequency  analysis  verified  it  as  a  TS.  Using  the  resulting  FCM  as  input, 
an  IRC  calculation  was  performed  in  the  direction  of  products.  The  walk  would  not  converge  due 
to  apparent  numerical  problems.  This  was  seen  by  the  fact  that  after  walking  successfully  for  six 
optimized  points  along  the  IRC,  the  bond  lengths  for  R(OH)  and  R(NN)  went  from  the  reasonable 
values  of  0.97  A  and  1.21  A,  respectively,  to  the  unrealistic  values  of  0.77  and  0.69  A,  respectively. 

The  calculations  were  repeated  for  the  determination  of  the  TS  structure  14a,  followed  by 
determination  of  the  FCM  for  the  optimized  critical  point.  Both  calculations  used  a  tighter 
convergence  criterion  on  the  SCF  density  matrix.  For  the  IRC  calculations,  Cartesian  coordinates 
were  used  to  follow  the  path  (instead  of  the  default  mass-weighted  internals),  and  the  convergence 
tolerance  on  the  SCF  density  matrix  was  set  to  “Converge =9”  from  the  standard  “Converge =8.” 
This  guarantees  an  RMS  error  in  the  density  matrix  of  1  x  10'9  or  less  and  an  absolute  error  in  the 
density  matrix  elements  of  1  x  10‘7  or  less.  The  default  settings  were  used  for  the  accuracies  of  all 
other  calculated  properties.  This  IRC  did  indeed  walk  successfully  to  a  converged  structure  (with 
default  gradient  convergence  criteria),  which  corresponded  to  MMLM’s  cis-planar  structure  7a  (our 
structure  K) — not  the  trans-planar  structure  7b  reported  by  MMLM.  There  were  no  indications  that 
any  gradient  step  produced  unrealistic  structural  parameters  as  in  the  previous  ran  using  the 
mass-weighted  internal  coordinates  and  default  SCF  convergence  criterion.  For  the  sake  of 
completeness,  the  QCISD/6-3 1  lG(d,p)  structure  corresponding  to  MMLM’s  7b  was  located  and  is 
shown  in  Figure  4  as  structure  P.  Both  geometries  are  planar,  and  the  most  notable  quantitative 


19 


difference  is  in  the  N-N-0  angle,  which  differs  by  10°  between  the  MP2(137°)  and  QCISD(127°) 
values.  These  results  seem  to  clearly  demonstrate  that  MMLM’s  cis-(H-O-N-N)  transition  state 
structure  14a  actually  connects  with  the  minimum  intermediate  7a  at  the  MP2  level,  which  is  also 
cis-(H-O-N-N)  planar  and  not  7b  as  previously  reported.  This  corrected  connection  of 
cis-(H-O-N-N)  transition  state  with  a  cis-(H-0-N-N)  minimum  structure  at  the  MP2  level  is  now 
consistent  with  the  QCISD/6-31G(d,p)  results. 

3.3.2  ly2-Hydrogen  Shifts.  It  is  conceivable  that  products  could  be  formed  through  a 
1 ,2-hydrogen  shift  from  trans-structure  G  on  Path  2.  This  section  will  describe  our  attempts  to  locate 
transition  states  corresponding  to  such  a  mechanism.  MMLM  report  two  1,2-hydrogen  shift  TSs 
corresponding  to  this  reaction  step:  4'  (nonplanar)  and  4  (planar).  These  structures  are  illustrated 
in  Figure  4  to  facilitate  comparison  with  the  current  results. 

The  nonplanar  4’  TS  was  sought  at  the  QCISD/6-3 1  G(d,p)  level.  The  structure  used  as  an  initial 
guess  corresponded  to  our  structure  C,  except  with  a  1 80°  rotation  about  the  N-N  bond.  Structure 
C  is  the  transition  state  for  the  1,2-hydrogen  shift  from  cis-structure  B.  An  SCF-level  FCM 
corresponding  to  the  initial  structure  was  used  as  an  initial  guess.  The  results  converged  to  a 
nonplanar  structure  denoted  as  N'  in  Figure  4.  Structure  N'  has  an  energy  of  60.38  kcal/mol  (see 
Table  3)  relative  to  reactants  and  was  verified  as  a  TS  through  vibrational  analysis  (see  Table  A-2). 
The  QCISD  structure  has  an  N-N  bond  length  that  is  0. 1 45  A  longer  than  that  of  structure  4'  reported 
by  MMLM.  Also,  the  QCISD  O-N-N  angles  are  considerably  smaller  than  those  of  structure  4' 
reported  by  MMLM.  The  QCISD  structure  has  an  05-N4-N2  angle  of  1 10°,  which  is  smaller  than 
the  MMLM  structure  4 '  value  by  22 0 .  The  QCISD/ 6-31  G(d,p)  structure  N '  more  closely  resembles 
the  planar  structure  4  reported  by  MMLM  except  for  the  nonplanarity  of  N'. 

These  structural  differences  could  be  due  to  differences  in  either  the  level  of  correlation  or  basis 
set.  To  test  the  basis  set  dependency,  the  transition  state  search  was  rerun  at  the  QCISD/6-3 1  lG(d,p) 
level,  and  the  resulting  converged  transition  state  structure  (confirmed  through  vibrational  analysis) 
was  almost  identical  to  the  structure  determined  using  the  smaller  basis  set  (structure  N').  The  bond 
angles,  bond  lengths,  and  dihedral  angles  of  the  QCISD/6-3 1 1  G(d,p)  structure  differed  by  1  %  or  less 
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from  N'.  Structure  N'  represents  the  only  geometry  located  in  this  study  associated  with  4',  with 
rather  poor  agreement  in  angles  between  the  two  structures  (see  Figure  4). 

MMLM’s  structure  4  was  reproduced  at  the  MP2/6-311G(d,p)  level  and  confirmed  as  a  TS 
through  vibrational  analysis.  A  subsequent  transition  state  search  for  structure  4  at  the 
QCISD/6-3 1  lG(d,p)  level  was  attempted  using  the  MP2/6-3 1 1  G(d,p)  optimized  structure  for  4  and 
a  QCISD/6-31  lG(d,p)  FCM.  The  converged  structure  is  a  TS  (see  Table  A-2)  and  is  denoted  as  N 
in  Figure  4.  It  is  markedly  different  from  MMLM’s  structure  4,  most  notably  in  the  N-N  bond 
separation.  The  N-N  bond  distance  in  N  is  1 .244  A,  which  is  considerably  shorter  than  the  1 .5 1 9  A 
of  structure  4.  The  ZPE  corrected  energy  of  N  relative  to  reactants  is  73.98  kcal/mol  (see  Table  3). 
And  once  again,  the  two  O-N-N  angles  for  structure  N  deviate  considerably  from  those  of  structure  4. 
Once  again  the  structural  dependence  on  basis  set  was  tested  by  searching  for  structure  4  at  the 
QCISD/6-3 1  G(d,f)  level  using  the  QCISD/6-3 1  lG(d,p)  FCM  and  geometry  as  the  initial  guess.  This 
calculation  converged  on  a  transition  state  structure  that  was  very  similar  to  N.  All 
QCISD/6-3  lG(d,p)  structural  parameters  agree  to  within  1%  or  less  of  the  QCISD/6-3  llG(d,p) 
values,  once  again  demonstrating  the  suitability  of  the  double-zeta  basis  set  for  structure 
determination. 

The  substantial  dissimilarities  between  the  MP2  and  QCISD  structures  for  4  vs.  N  and  4'  vs.  N' 
raise  suspicions  about  the  suitability  of  MP2  for  predicting  critical  points  on  the  PES  for  this  system. 
From  the  standpoint  of  energetics,  the  QCISD/6-3 1  lG(d,p)  barriers  corresponding  to  the  planar  and 
nonplanar  structures  (N  and  N',  respectively)  are  74  and  60  kcal/mol  (Table  3)  relative  to  reactants. 
Since  these  were  so  much  higher  in  energy  than  the  rate-limiting  step  on  Path  2,  no  further  studies 
were  conducted  in  search  of  a  reaction  path  through  the  trans-1,2  hydrogen  shift. 

One  structure  on  MMLM’s  MERP  that  has  not  yet  been  addressed  by  this  study  at  the  QCISD 
level  is  structure  7b.  It  has  been  verified  through  an  IRC  walk  that  TS  N  connects  to  a  minimum 
structure  qualitatively  similar  to  7b  at  the  QCISD/6-3  llG(d,p)  level.  The  IRC  from  N  completed 
17  steps,  at  which  point,  an  optimization  was  run  that  produced  the  minimum  structure  P  in  Figure  4. 
The  relative  energy  of  P  is  given  in  Table  3  as  6.02  kcal/mol,  nearly  identical  to  the  relative  energy 
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of  6.1  kcal/mol  for  7b  reported  by  MMLM.  This  route  to  forming  P  (MMLM’s  7b)  via  TS  N  does 
not  preclude  the  possibility,  in  fact  the  likelihood,  of  forming  P  through  the  path  predicted  by 
MMLM  (i.e.,  TS  14a  -  7b  [our  P],  which  occurs  by  a  rotation  about  the  N-O  bond  in  N-O-H 
[Figure  3]).  A  comparison  of  P  with  7b  shows  that  the  most  noteworthy  structural  difference 
between  the  QCISD  and  MP2  geometries  is  a  10°  difference  in  the  angle  N2-N4-05. 

4.  Summary  and  Conclusions 

We  have  presented  results  from  an  ab  initio  study  of  the  reaction  HNO  +  NO  -  N20  +  OH. 
Geometries  of  critical  points  were  optimized  at  the  QCISD/6-31G(d,p)  level  and  characterized 
through  normal  mode  analyses.  The  results  of  this  study  provide  two  paths  leading  to  formation  of 
products  through  formation  of  various  reaction  intermediates.  The  reaction  paths  were  established 
through  IRC  calculations  leading  from  the  transition  states  to  connecting  minima,  which  resulted  in 
what  appears  to  be  the  MERP  and  a  second,  higher-energy  pathway.  No  claim  is  being  made  that 
these  are  the  only  pathways  leading  to  products  of  Reaction  I.  Rather,  they  constitute  two  completely 
connected  pathways.  The  transition  state  structure  corresponding  to  the  highest  energy  barrier  on 
the  MERP  (denoted  as  H)  in  this  study  supports  the  conclusion  of  MMLM  that  their  barrier 
structure  3  represents  the  rate-determining  step  in  the  reaction  HNO  +  NO  -  N20  +  OH. 
Establishment  of  this  barrier  height  by  MMLM  is  a  very  important  result,  since  it  is  the  first  strong 
piece  of  evidence  supporting  the  kinetic  modelers  ’  hypothesis  that  this  reaction  can  play  a  major  role 
in  the  gas  phase  chemistry  of  burning  nitramine  propellants.  MMLM  also  used  QCISD,  but  only  in 
determining  the  geometries  of  the  reactants,  products,  and  TS  3.  Our  best  ZPE  corrected  energy  for 
this  barrier  height  is  24  kcal/mol,  to  be  compared  with  MMLM’s  “best”  value  of  21 .6  kcal/mol  and 
the  value  26  ±  5  kcal/mol  extracted  from  experimental  data. 

Apart  from  agreeing  on  the  rate-determining  barrier,  the  current  results  differ  most  prominently 
from  those  of  MMLM  in  the  definition  of  the  remaining  MERP.  hi  particular,  our  MERP  entrance 
barrier  has  a  trans-(O-N-N-O)  structure  rather  than  the  cis-structure  reported  at  the  MP2  level  by 
MMLM.  Reaction  proceeds  through  formation  of  a  trans-minimum  intermediate  G,  which  then  must 
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overcome  the  rate-determining  barrier  H.  Transition  state  structure  H  is  trans-,  in  qualitative 
agreement  with  MMLM’s  MP2  structure  3.  After  crossing  barrier  H,  the  reaction  proceeds  to 
products,  while  retaining  the  relative  cis-(H-0-N-N)  and  trans-(O-N-N-O)  torsional  angle 
relationships  seen  in  H.  In  contrast,  MMLM’s  TS  cis-14a  (our  J)  proceeds  to  a  minimum 
intermediate,  which  is  trans-(H-O-N-N),  and  then  over  their  exit  barrier  13,  which  is  still 
trans-(H-O-N-N)  to  products.  While  their  exit  barrier  13  is  lower  than  our  cis-(H-O-N-N)  barrier, 
the  IRC  from  our  J  (MMLM’s  14a)  leads  to  a  cis-structure  at  both  MP2  and  QCISD — not  the  trans¬ 
structure  in  MMLM.  Almost  certainly,  there  exists  a  transition  state  (probably  low  energy)  leading 
from  either  I  or  K  to  such  a  trans-(H-O-N-N)  minimum  intermediate  like  their  13;  however,  it 
appears  that  14a  is  not  the  correct  candidate.  The  current  study  predicts  an  MERP  connection  of 

Reactants -F-G-H-I-J-K-L-  Products, 
whereas  MMLM  predicts  (using  our  nomenclature  for  structures) 

Reactants  -*A-»G->H->I-J-,P-'E-*  Products. 

Some  of  the  QCISD  structures  are  quantitatively  and  qualitatively  dissimilar  to  the  MP2  results 
of  MMLM  and  Mebel,  Morokuma,  and  Lin  [7].  For  example,  the  planar  transition  state  6'  (for 
rotation  about  the  N-N  bond)  reported  by  Mebel,  Morokuma,  and  Lin  [7]  exists  not  as  a  TS  but  as 
a  minimum  structure  at  the  QCISD/6-31  lG(d,p)  level.  Its  QCISD  analogue  in  this  study  is  the  TS 
structure  M,  which  is  decidedly  nonplanar.  Based  on  its  energy  relative  to  reactants,  M  could  most 
certainly  be  involved  in  the  reaction  dynamics  that  lead  to  products.  Due  to  the  expense  of  the 
QCISD/6-3 1  lG(d,p)  calculations,  only  a  few  other  such  comparisons  for  structures  not  on  Path  1  or 
Path  2  were  pursued,  including  the  trans-(O-N-N-O)  1,2  H-shift  TS  structures  4  and  4'  of  MMLM. 
The  QCISD  did  indeed  find  analogues  to  these  transition  states,  labeled  N  and  N'  in  Figure  4.  But 
the  R(N-N)  bond  lengths  differ  between  4  and  N  by  0.28  A  and  one  O-N-N  angle  by  20°.  Similar 
disparities  exist  between  the  4'  and  N'  geometries,  where  the  R(N-N)  differs  by  0.14  A  and  one 
O-N-N  angle  by  21  °.  Fortunately,  under  combustion  conditions,  N'  and  N  are  not  likely  to  play  a 
role  in  the  reaction  leading  to  products.  Neither  MMLM  nor  Mebel,  Morokuma,  and  Lin  [7]  report 
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a  transition  state  structure  analogous  to  our  exit  barrier  L,  a  structure  with  cis-(HONN)  and 
trans-(ONNO).  One  last  notable  difference  between  our  QCISD/6-31G(d,p)  calculations  and  the 
MP2  results  reported  by  Mebel,  Morokuma,  and  Lin  [7]  is  for  our  cis-planar  structure  B.  Mebel, 
Morokuma,  and  Lin  [7]  report  a  cis-planar  structure,  labeled  6',  which  is  similar  to  our  minimum 
structure  B,  but  they  report  that  it  corresponds  to  the  transition  state  for  the  internal  rotation 
around  the  N-N  bond.”  However,  our  structure  B  is  a  minimum,  not  a  transition  state.  These 
predicted  differences  in  structures  between  QCISD  and  MP2,  some  of  which  are  quite  large,  could 
produce  subsequently  large  errors  in  the  relative  energetics,  vibrational  frequencies,  and  normal 
modes  of  the  species  on  the  PES.  This,  coupled  with  the  differences  in  connectivity  predicted  by  the 
IRC  walks  could  play  a  major  role  in  the  outcome  of  a  dynamical  simulation  for  Reaction  I. 
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Table  A-l.  QCISD/6-31G(d,p)  Vibrational  Frequencies  (in  cm'1)  for  Critical  Points  on  the 
PES 


HNO 

NO 

N20 

OH 

A 

B 

C 

1580 

1950 

590 

3778 

595/ 

284 

2209/ 

1637 

590 

192 

351 

238 

2982 

1292 

259 

490 

327 

2342 

569 

756 

391 

710 

926 

492 

1411 

1381 

850 

1446 

1474 

1170 

1982 

1602 

1715 

3191 

3557 

2203 

C' 

D 

E 

F 

G 

H 

I 

2004 / 

285 

787/ 

659/ 

231 

654/ 

229 

107 

394 

208 

127 

418 

248 

419 

308 

562 

213 

241 

629 

603 

654 

342 

773 

552 

476 

629 

801 

663 

533 

927 

639 

806 

1033 

1024 

879 

834 

1155 

1083 

1401 

1379 

1191 

1192 

1187 

1416 

1272 

1414 

1494 

1280 

1410 

1773 

1627 

1947 

2002 

1569 

1566 

1513 

2773 

3878 

3831 

3112 

3427 

2108 

3622 

J 

K 

L 

M 

N 

N' 

2111/ 

341 

837/ 

212/ 

2290/ 

1888/ 

290 

378 

42 

358 

161 

179 

448 

403 

281 

462 

383 

372 

577 

641 

435 

674 

468 

381 

751 

1029 

589 

862 

520 

591 

1141 

1228 

1142 

1334 

871 

1117 

1233 

1461 
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3658 
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Table  A-2.  QCISD/6-311G(d,p)  Vibrational  Frequencies  (in  cm'1)  for  Critical  Points  on  the 
PES 
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